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We propose a concept of superconductivity controlled by orbital degree of freedom taking 
CeMIns (M — Co, Rh, and Ir) as typical examples. A microscopic multiorbital model for 
CeMIns is analyzed by fluctuation exchange approximation. Even though the Fermi-surface 
structure is unchanged, the ground state is found to change significantly among paramagnetic, 
antiferromagnetic, and d-wave superconducting phases, depending on the dominant orbital 
component in the band near the Fermi energy. We show that our picture naturally explains the 
different low-temperature properties of CeMIns by carefully analyzing the crystalline electric 
field states. 

KEYWORDS: crystalline electric field, heavy-term ion superconductor, fluctuation exchange approximation, 
thermal expansion 



Since the discovery of heavy-fermion superconductiv- 
ity in CeCu2Si2, 1 it has been one of the central issues 
in the research field of condensed matter physics to un- 
veil unconventional superconductivity in strongly corre- 
lated electron systems. In particular, it is important to 
determine a key parameter for controlling the appear- 
ance of superconductivity. Among heavy-fermion super- 
conductors, CeMIns (M = Co, Rh, and Ir) have po- 
tential for systematic understanding of the mechanism 
and control parameter of superconductivity, since vari- 
ous ordered phases have been found in the same crystal 
structure. At ambient pressure, a relatively high super- 
conducting transition temperature T c = 2.3 K has been 
reported for CeCoIns, 2 whereas T c = 0.4 K for Celrlns. 3 
CeRhIn5 is an antiferromagnet with a Neel temperature 
Tn = 3.8 K at ambient pressure, but it becomes super- 
conducting with T c ~ 2 K at pressure p > 1.6 GPa. 4 

The observed power-law temperature dependences of 
specific heat, 2,5 ' 6 thermal conductivity, 5 and spin-lattice 
relaxation rate 7-10 suggest the presence of line nodes in 
the superconducting gap function. In addition, the four- 
fold magnetic-field angular-dependences of thermal con- 
ductivity 11 and specific heat 12 clearly indicate c£-wave 
superconductivity for CeCoIns. The existence of anti- 
ferromagnetism in CeRhlns at ambient pressure is also 
consistent with d-wave superconductivity induced by an- 
tiferromagnetic fluctuations. 

Concerning the control parameter of superconductiv- 
ity, one may first consider carrier density, as in the case 
of high-T c cuprates. In addition, Fermi-surface topology 
also plays a crucial role in the occurrence of supercon- 
ductivity and the determination of the symmetry of the 
gap function. However, the carrier densities are essen- 
tially the same among CeMIns, because /-electron num- 
ber per Ce ion is not changed by M. Furthermore, Fermi 
surfaces observed by de Haas-van Alphen experiments 
and band-structure calculations are similar among these 
compounds. 13-18 Thus, neither the carrier density nor 
Fermi-surface topology is the main control parameter of 
superconductivity and antiferromagnetism in CeMIns. 

Here, we point out another important ingredient, the 



crystalline electric field (CEF) effect, in CeMIns. Imme- 
diately after the discovery of CeMIns, Takimoto et al. 
have stressed the importance of the CEF effect in con- 
trolling superconductivity, 19 but only level splitting was 
taken into account. Note that CEF potential affects not 
only level splitting, but also the CEF ground-state wave- 
function. In fact, recent inelastic neutron scattering ex- 
periments have revealed that level splittings are almost 
the same among CeMIns, whereas CEF wave-functions 
are quite different. 20 The importance of orbital states has 
been discussed for Na x Co02-yH 2 0, 21 but the effects of 
the change in orbital states on superconductivity have 
not been studied systematically so far. It is conceptu- 
ally important to clarify superconductivity controlled by 
orbital states. 

In this Letter, we apply fluctuation exchange (FLEX) 
approximation to an /-electron model constructed on 
a square lattice for CeMIns. Such perturbative the- 
ories have been applied to one- /-orbital models for 
CeMIns. 22-24 In this study, we focus on the effect of a 
CEF using a model including all the states with the to- 
tal angular momentum j — 5/2, which are split into one 
Tq and two V-? doublets under a tetragonal CEF. Among 
them, the wave functions of T7 states sensitively depend 
on CEF potential. We change the wave functions fixing 
level splittings, and determine both T c and Tn within 
FLEX approximation. Even though the Fermi-surface 
structure does not change, the ground state changes de- 
pending on the T7 wave-functions among the param- 
agnetic, antiferromagnetic, and d-wave superconducting 
states. We show that the obtained phase diagram is rel- 
evant to CeMIn.5. 

Using the /-electron basis under a cubic CEF for con- 
venience, we consider a three-orbital model based on a 
j-j coupling scheme on a square lattice given by 
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tron at site r with the pseudospin a and the orbital r, 
n rTa = c\. rr7 c rTa , and n rT — \\^ a n rT<J . Note that a (=\ 
and |) distinguishes two states in each Kramers doublet. 
On the other hand, r is introduced to distinguish three 
kinds of Kramers doublets under a cubic CEF. Here, a 
and P denote two r§, while 7 indicates r 7 . We use t^ T , 
for the hopping amplitude between the r' state at site 
r + fi and the r state at r, where n is a vector connecting 
nearest-neighbor sites. We consider the hopping through 
(ffcr) bonding, given by = iSf' = -VStffl = 

-V3t£ 0) = V3t% a) = V3tf a a) = 3^ 0) = 3$ o) = t 
and zero for the other cases, 25 where t = 9(//er)/28 and 
a is the lattice constant. The coupling constants U and U' 
denote intra- and inter-orbital Coulomb interactions, re- 
spectively. Note that we ignore the Hund's rule coupling, 
pair-hopping interaction, and the other interactions for 
simplicity. Then, owing to the symmetry requirement, we 
should set U — U' . We also note that t,U, and V should 
be regarded as renormalized ones for heavy electrons. 26 
The tetragonal CEF term -ffcEF is given by 

Hqef = Y.^°t + BlOt + BjOl), (2) 

r 

where B^ is the CEF parameter and 0^ r is the 
Stevens operator equivalent at site r. 27 We can rewrite 
cq. (2) with the basis of the cubic CEF eigenstates as 
#cef = Z)r,<7,T,T' B TT'ct TtT c rT > a , where B TT > is given 
by an appropriate linear combination of B^. From the 
diagonalization of Hqef, we express the CEF parame- 
ters as B® = [-4A 67 + A 7 (2cos(2<9) - v/5sin(20))]/84, 
B\ = [3A 67 + A 7 (2cos(26») - V5sin(26»))]/1260, and 
B\ = V5A 7 (V5cos(26l) + 2sin(26»))/360, where A 67 and 
A 7 determine CEF energy levels, while 9 character- 
izes eigenstates under a CEF potential. Then, the CEF 
energy levels are — A@ 7 /3 — A 7 /2 for the states, 
(4 7a cos0 + 4 a(7 sin0)|O), -A 67 /3 + A 7 /2 for the T 7 2) 
states, (— cJ. 7<T sm8 + c) raa cos0)|O), and 2Ag 7 /3 for the 

Tq states, cj^jO), where |0) denotes vacuum. 

Now we apply FLEX approximation, which has been 
extended for multiorbital models. 28-31 Since a is a con- 
served quantity in the present model, the Green's func- 
tion does not depend on a and is represented by a 3 x 3 
matrix. In a matrix form, the Dyson-Gorkov equation 
is given by G{k) = G(°>(fc) + G(°> (jfe)E(jfe)G(fc), where 
G(k) is the Green's function and G^°\k) is the nonin- 
teracting Green's function. Here, we have introduced the 
abbreviation k = (fc, ie„), where k is the momentum and 
e n = (2n + l)nT is the fermion Matsubara frequency 
with an integer n and a temperature T. The self-energy 
is given by 

Z TlT2 (k) = Jj^2 V riri:r 2 Ti(q)G TiTi (k- q), (3) 

where N is the number of lattice sites, q = (q,iuj m ), 
and u) m = 2mirT is the boson Matsubara frequency. 
The fluctuation exchange interaction is given by V{q) = 
l[U s X s {q)U s - U s X [a) {q)U s /2 + U s ] + \[U c X c {q)U c - 
U c x {0) (qW c /2 - U c ). The matrices U s and U c are given 
bv II s = U c = U II s , , = -U c , , = U' 

U J t t \ t t u tt;tt u i u tt';tt' ° TT> ;tt' u ' 



and U^ T . T , T , — 2U' , where r ^ r', and the other matrix 
elements are zero. The spin and charge parts of the sus- 
ceptibility are given by x s (<?) = X^°Hq)[^ ~ U s x (QS> (q)] -1 
and X c (q) — x'°H<z)[l + U C X ^ (q)] -1 , respectively, where 
XTiT 2 ;T 3 T4,(q) — — jj X^fe G TlT3 (k + q)G TiT2 (k). 

In FLEX approximation without vertex corrections, 
the magnetic susceptibility is given by 32 

Xv{q) X ^ ' XtiT2;T3T4 iq)^T2CT2;T 1 (J 1 ^ T 3 <J 3 \T&a& ® <7\d1 ®<7i<7 3 ' 

TI-T4 
cr 1 -CT4 

(4) 

where v — x, y, or z, cr are the Pauli matrices, and 
J^ 2(T2 . Tl(Tl is the matrix element of the operator of the 
dipole moment. In this paper, we renormalize J v so that 
the sum of squares of eigenvalues is unity, for conve- 
nience. The linearized equation for the anomalous self- 
energy <f> is expressed as 

= -5 E (5) 

where ^ denotes the pseudospin singlet (S) or triplet (T) 
pairing state and the kernel is given 

Er^VAr^r^ ~ k')G T:3T[ (k')G T ^(-V) . The effec- 
tive pairing interactions are given by V s (q) = 
\[U s x s {q)U s + U S /2]-\[U c x c {q)U c -U c /2] and V T (q) = 
- \ [U s x s (q) U s + U s /2]-\ [U c x c (q) U c - U c /2] . The tran- 
sition temperature of superconductivity is given by the 
temperature where eq. (5) has a nontrivial solution. 

In the following, we show results for a 32 x 32 lattice 
for U = U' = 4t, A 67 = 1M, and A 7 = t. The results are 
not so sensitive to precise values of Ag 7 and A 7 as long 
as they are in the order of t. In the calculation, we use 
1024 Matsubara frequencies. The number of / electrons 
per site is fixed to be one, since we are considering Ce 3+ 
ions. The present model with the CEF parameter 6 is 
invariant for 6 — > 6 + it. In addition, when U = U', the 
model is also invariant for 9 — ► —9. Thus, it is enough to 
consider < 9 < ir/2. 

Figures l(a)-l(c) show the magnetic susceptibility 
Xv(q) = Xv(q,iu m = 0) for 6 = 0, 7r/4, and tt/2, re- 
spectively. It is observed that Xv(q) strongly depends on 
9. In the following, we explain this dependence from the 
viewpoint of itinerant and localized orbitals. 

First, let us discuss the noninteracting susceptibility 
X^\q) [Figs. l(d)-l(f)]. In Figs. l(g)-l(l), we show non- 
interacting band structures and Fermi-surface curves. 
Even if 9 is changed, the band that crosses the Fermi 
level E-p has similar dispersion at Ep, leading to an al- 
most the same Fermi surface. Note, however, that in gen- 
eral, Xv*\q) depends on the band structure and is not 
determined solely by the Fermi surface structure {e.g., 
nesting property), particularly for multiorbital systems. 
The nearly flat band composed mainly of the 7 orbital 
locates near Ep for small 8, 33 suggesting that electrons 
near Ep have nearly localized character. 34 Thus, X ^ (?) 
for 9 — is almost flat in momentum space. For small 9, 
the magnitude of X ^ (q) is large, since x^ ~ Pa at suffi- 
ciently low T, where po is the density of states at Ep, and 
it becomes large when the localized character is strong. 
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Fig. 2. Phase diagram for U = U' = At, A 6 7 = 1.5t, and A 7 = t. 
Solid squares denote antiferromagnetic (AF) transition tempera- 
tures, while solid circles denote superconducting (SC) transition 
temperatures with Bi g symmetry. 
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Fig. 1. (a)-(c) Magnetic susceptibility at U = U' = At and T = 
0.02t, for = 0, 7r/4, and 7r/2, respectively, (d)-(f) Magnetic 
susceptibility at U = U' = and T = 0.02t, for = 0, tt/4, and 
7r/2, respectively, (g)— (i) Energy band structures at U = U' = 
for 9 = 0, 7r/4, and tt/2, respectively. The Fermi energy is set 
at zero (thin lines), (j)— (1) Fermi-surface curves at U = U' = 
for 6 = 0, 7r/4, and n/2, respectively. The parameters are set at 
A 67 = 1.5t and A 7 = t. 
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With increasing 6, the localized character is weakened. 
Thus, Xv'Hl) shows a structure in q space and becomes 
smaller in magnitude. In particular, we observe a moder- 
ate enhancement around q = Q = (it / 'a, it /a) for 9 = it /A 
and tt/2. 

Now, we consider the effect of on-site Coulomb interac- 
tions. We expect that the Fermi surface is not markedly 
changed by the Coulomb interaction. First, note that the 
on-site Coulomb interactions suppress the qr-independent 
part of (<?) due to the correlation effect beyond ran- 
dom phase approximation (RPA). Thus, for small 9, 
Xv(q) is totally suppressed in comparison with that of 
RPA. For larger 9, the structure in X^(q) grows to a 
peak at q = Q in Xv{q)- On the other hand, such a peak 
becomes high and sharp, when the Coulomb interactions 
are effectively strong as for almost localized electrons at 
small 9. Thus, due to the combination of the band struc- 
ture and the effect of the Coulomb interactions, XvW) 
is large for a moderate value of 9 ~ tt/4 with a peak at 

q = Q 

In Fig. 2, we show a T-9 phase diagram. In this paper, 
we define Tn as a temperature at which Xz{Q) reaches 
100/t. Note that the susceptibility does not become as 
large as 100/i for the other components Xx and Xyi 01 
for q other than Q. The overall structure in Fig. 2 is con- 
sistent with the 9 dependence of the magnetic property. 
Note that a c?-wave singlet superconducting state with 
Bi 9 symmetry appears at 9 > 0.297T next to the antifer- 
romagnetic phase. On the other side of the antiferromag- 
netic phase, T c is very low, 35 even if the superconducting 
phase exists. This is consistent with the fact that the en- 
ergy scale in the small- 9 region becomes small owing to 
the strong localized character of electrons. 



Fig. 3. Temperature dependences of (a) a®, (b) aJJ, and (c) — a\ 
for single ion under CEF potential deduced from neutron scat- 
tering experiment on CcRhlns. 20 Solid curves denote the values 
without a magnetic field and dashed curves denote those in a 
magnetic field H z = 18 T along the c-axis. 



Now we discuss possible relevance of our results to 
CeMIns. For this purpose, it is necessary to estimate 
9 for each material. Among the CEF parameters, the 
sign of B\ is quite important, since the orbital state is 
markedly affected by the sign. Although inelastic neutron 
scattering is a powerful method of determining CEF en- 
ergy levels, the sign of B\ cannot be determined solely 
by neutron scattering experiment. Thus, we should per- 
form CEF analysis of a quantity sensitive to (Of), where 
(• • • ) denotes the expectation value in terms of H CEF . 
Since the mode of charge distribution corresponding to 
finite (Of) couples to the lattice, thermal expansion is 
a good quantity for determining the sign of B\. Ther- 
mal expansion has been measured for CeMIns, but the 
analysis focused only on the second-order term (O!,)- 36-40 
It is necessary to include fourth-order terms to analyze 
thermal expansion, as in the case of cubic systems. 41 

Figures 3(a)— 3(c) show the temperature dependence of 
a° 2 = d(0°)/dT, a° = d(Ol)/dT, and a\ = d(Oj)/dT, 
respectively, for the CEF parameters deduced from the 
neutron scattering experiment for CeRhlns, 20 assuming 
B\ > 0. A thermal expansion coefficient is given by a 
linear combination of a®, ot% and a\. For B\ < 0, a\ 
changes its sign, whereas a° an d ce® do not. Although 
we do not show the calculated results for CeCoIns and 
Celrlns, it has been found that the overall features are 
similar among CeMIns. Figures 3(a)-3(c) show that the 
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Fig. 4. Schematic phase diagram for CeMIns. Solid square and 
circles denote and T c , respectively. Curves are guides for the 
eyes. See the main text for the values of 6. 



magnetic-field effect on aP> and a% is strong, but weak 
on a\. From the experimental results, 38,39 the magnetic- 
field effect on the thermal expansion coefficient a c along 
the c-axis for CeRhIn 5 is not so significant as that ob- 
served in Figs. 3(a) and 3(b). Rather it looks similar to 
that in Fig. 3(c). Moreover, the minimum with a negative 
value in —a\ shifts toward higher temperatures under a 
magnetic field, as is experimentally observed in a c of 
CeRhIn5. Thus, we highly expect that a c in CeMIns is 
mainly determined by a\. For CeCoIns and Celrln5, a c 
is always positive at the temperatures reported, 36,37 in 
sharp contrast to that for CcRhlns. These observations 
indicate that the signs of B\ for CeCoIn 5 and Cclrins 
are the same, but are different from that for CeRhIn 5 . 

For further quantitative analysis, it is necessary to 
have more precise knowledge on elastic constants and/or 
magnetoelastic coupling, but it is out of the scope of this 
paper. Rather, we phenomcnologically assume B\ > for 
CcRhIn 5 . Then, B\ is negative for CeCoIn 5 and Celrln 5 . 
By using the CEF parameters deduced from the neu- 
tron scattering experiment, 20 we obtain 9 = 0.16-7T for 
CeRhIn 5 , 9 = -0.30tt for CeCoIn 5 , and 9 = -0.38tt for 
Celrlns. In Fig. 4, we show a schematic phase diagram 
using these values of 9 and experimental results for T c 
and Tn. This is consistent with the theoretically deter- 
mined phase diagram in Fig. 2. In addition, CcRhIn 5 
locates in a small-# region, indicating that / electrons 
in CeRhIn 5 have a relatively localized character, as is 
suggested from de Haas- van Alphen experiment. 17 Thus, 
we conclude that the characteristics of CeMIn5 can be 
captured by our model and that the main control pa- 
rameter of antiferromagnetism and superconductivity for 
CeMIns is 9 which describes CEF wave-functions. To 
confirm this, it is important to experimentally determine 
CEF parameters for CeMIn 5 under pressure. 

In summary, we have studied the /-electron model in- 
cluding all the states with j = 5/2 by FLEX approx- 
imation. We have found three kinds of ground state, 
i.e., paramagnetic, antifcrromagnetic, and d-wave super- 
conducting states, by changing the CEF wave- functions 
characterized by the parameter 9, even if we fix level 
splitting. Such phase change originating from the differ- 
ence in character of /-electron orbitals explains well the 



difference in CeMIn5. This finding shows that, in general, 
orbital character can be a control parameter of super- 
conductivity, in addition to Fermi surface topology and 
carrier density. The orbital-controlled superconductivity 
would be realized also in other /-electron materials with 
active orbital degrees of freedom, and the present study 
will provide an important step toward the investigation 
of such exotic superconductivity. 
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